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Abstract. We propose an effective algorithm that enumerates (and 
actually finds) all 3-edge colorings and Hamiltonian cycles in a cubic 
graph. The idea is to make a preliminary run that separates the vertices 
into two types: "rigid" (such that the edges incident to them admit a 
unique coloring) and "soft" ones (such that the edges incident to them 
admit two distinct colorings), and then to perform the coloring. The 
computational complexity of this algorithm is on a par with (or even 
below) the fastest known algorithms that find a single 3-edge coloring 
or a Hamiltonian cycle for a cubic graph. 



1. Preliminaries 

Let us recall here some basic facts about the relationship between 3- 
edge colorings (also called Tait colorings) and Hamiltonian cycles in cubic 
graphs; the details can be found in pQ. By a cubic graph we understand a 
connected 3-regular multi-graph that is allowed to have double edges, but 
no loops (obviously, a graph with loops cannot have either edge colorings or 
Hamiltonian cycles) . Consider a set of three distinct elements called "colors" 
(say, {r, g, b}, where r stands for "red" , g — for "green" , and b — for "blue"). 
A 3-edge coloring, or Tait coloring is an assignment of a color to every edge 
such that the edges incident with each vertex have distinct colors 1 . Every 
3-edge coloring of a cubic graph G gives rise to three distinct 2-factors (that 
is, 2-regular spanning subgraphs) of G called Tait cycles: each Tait cycle is 
the union of edges painted in two colors out of the three (the complement to 
a Tait cycle is a perfect matching - the union of disjoint edges painted in the 
third color). If a cubic graph G has a Hamiltonian cycle, then G also admits 
a 3-edge coloring that is unique up to a permutation of colors: just paint 
the Hamiltonian cycle (which always has even length) in two intermittent 
colors, and paint the complement perfect matching in the remaining third 
color. 
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The above connection between 3-edge colorings and Hamiltonian cycles 
suggests the following method of enumerating (and actually finding) all the 
Hamiltonian cycles in a cubic graph: 

(1) Find all 3-edge colorings (Tait colorings) of a given cubic graph up 
to permutations of colors; 

(2) Find all the corresponding 2-factors (Tait cycles); 

(3) Check for connected 2-factors (Hamiltonian cycles). 

This procedure gives the complete list of Hamiltonian cycles in a cubic graph. 

For the enumeration of 3-edge colorings we use an exhaustive backtrack- 
ing algorithm that works in two runs. During the first run it dynamically 
separates all the vertices into two types: "rigid" ones (that admit a unique 
coloring of edges incident with them), and "soft" ones (with exactly two 
possibilities of coloring the incident edges); no backtracking is needed on 
this stage. The second run is the actual painting of edges: after it success- 
fully colored the graph or was unable to complete the coloring, it returns to 
the last visited soft vertex and tries a different possibility. The details are 
explained in the next section. 

2. Description of the algorithm 

First we partition the set of vertices V(G), of the graph G into two disjoint 
subsets V(G) = R(G)US(G) of rigid (R) and soft (S) vertices. We note that 
this partition is not canonical. Initially we put S = R = and dynamically 
change their content. We also introduce a temporary set U of unidentified 
vertices that we already visited, and an ordered list of colored edges E. 
We label the vertices of G by integers {0, . . . , n — 1}, n = #V(G). For 
the vertex with number i we denote the numbers of adjacent vertices by 
n-o < n \ < n 2- An edge connecting i and j we denote by We start 

at the vertex and add it to the set R of rigid vertices. We add the three 
edges [0, no], [0, n°], [0, n°] incident with it to the list E, and we add their 
endpoints n^n^n® to U. Now we check if any of the vertices in U are the 
endpoints of at least two edges in E. If this is the case, we move all such 
vertices from U to R, and for every such vertex we also add the remaining 
third edge incident to it to the set E (if it is not already there) . We continue 
the above procedure until there is no vertex left in U that is an endpoint 
of at least two edges in E. Now pick the vertex from U with the smallest 
number, say i, and move it into S. Note that i is an endpoint of a single 
edge in the current set E. Next, we append to E the two remaining edges 
incident with i, add their endpoints to U, and again check if any of the 
vertices in U bound at least two edges in E. If they do, such vertices are 
moved to R, and the missing edges incident with these vertices are added 
to E. Otherwise, we pick a vertex from U with the smallest number, move 
it to S and repeat the procedure until U becomes empty, or, equivalently, 
until E coincides with E(G). Since G is connected, this would mean that 
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Figure 1 



1-skeleton of a 3-cube. 



V(G) = R U S. Setting R(G) = R and S(G) = S we obtain the required 
partition. 

The final list of edges E provides the order in which we attempt to paint 
the edges of the graph. As above, we start at the vertex and paint the 
edges [0, no], [0, n°], [0,n[>] incident with it in colors a,b,c respectively. If 
the next edge in E has a rigid vertex as its endpoint, then there is a unique 
color left for the remaining third edge incident with that vertex. If the 
next pair of edges is incident with a soft vertex, we choose one of the two 
options for painting these edges (in case we visit this vertex for the first 
time). Continuing this way we may successfully reach the end of the list 
E and get a complete edge coloring that we save for our record. It may as 
well happen that the procedure ends prematurely when two edges incident 
to the same vertex are painted in the same color. In both cases we return 
in E to the previous soft vertex that was visited only once, and start over 
painting edges in a different way. At the end we get the list of all possible 
Tait colorings and check which of them produce connected Tait cycles (that 
is, Hamiltonian cycles). 

Let us now illustrate this algorithm on a simple example. 

Consider the graph G on 8 vertices {0,1,2,3,4,5,6,7} shown on Fig. 
1. We start with S = R = U = E = %. The initial vertex is a rigid 
vertex, R = {0}. Add the edges [0, 1], [0, 3], [0, 7] to E and the vertices 
1, 3, 7 to U. None of these vertices is an endpoint of two edges in E, so we 
remove the vertex 1 from U, append it to S, and add the edges [1,2] and 
[1,6] to E. Again, none of the vertices in U is an endpoint of two edges 
in E, so we move the vertex 2 to S and add the edges [2, 3] and [2, 5] to 
E. Now the vertex 3 bounds two edges in E, namely, [0,3] and [2,3]. We 
move it to R, add the edge [3, 4] to E and the vertex 4 to R. None of the 
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vertices {4, 5, 6, 7} in U bounds at least two edges, so we move the vertex 

4 to S and add the edges [4, 5] and [4, 7] to E. We see that the vertices 

5 and 7 become rigid, and so does the vertex 6. Thus, we get S{G) = 
{1,2,4}, R(G) = {0,4,5,6,7}, and the ordered set of edges is E(G) = 
{ [0, 1] , [0, 3] , [0, 7] , [1, 2] , [1, 6] , [2, 3] , [2, 5] , [3, 4] , [4, 5] , [4, 7] , [5, 6] , [6, 7] } . 

Let us start coloring the graph. We paint the edges [0, 1], [0, 3], [0, 7] in 
colors r, g, b respectively. There are two possible ways of coloring the edges 
incident with the first soft vertex 1. We first paint [1, 2] in color g and [1, 6] 
- in color b. At the next soft vertex 2 we again have two options, and 
choose the first one of them — paint [2, 3] in color r, and [2,5] — in color b. 
The edge [3, 4] incident with the rigid vertex 3 necessarily have color b. At 
the last soft vertex 4 both possibilities of coloring the edges [4, 5], [4, 7] lead 
to complete edge colorings of G. Now we return to the previous soft point 2 
and paint [2, 3] in b, and [2,5] — in r. This gives us one more Tait coloring 
(one of the options cannot be completed) . Finally, returning to the first soft 
point 1 and painting the edges [1,2], [1,6] in colors b,g respectively, we get 
the last coloring. Thus, there exist 4 distinct Tait colorings of the graph G 
(up to permutations of colors) listed in the following table: 
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The first one of these 4 edge colorings has no Hamiltonian cycles associ- 
ated with it, whereas the other three ones produce two Hamiltonian cycles 
each. They are listed below (we indicate in brackets the corresponding al- 
ternating colors): 
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3. Computational complexity 

It is clear that the computational complexity of this algorithm is of order 
2# S ( G ) . Since the set of soft vertices S (G) depends on the ordering of vertices 
of G, the complexity also depends on this ordering. To give an upper bound 
for #S(G) for a simple graph without double edges we note that, every time 
we add two new edges to E incident with a soft vertex, we encounter one of 
the three possibilities: 

(1) Both new endpoints belong to U; 

(2) One new endpoint belongs to U an one to R (the vertex from R then 
gives rise to new vertices that are added to U U R) ; 

(3) Both new endpoints belong to R. 

In any case, every soft vertex gives rise to at least two new vertices in U U R 
and the lower bound 2 is attained when the both endpoints of the edges 
incident with a soft vertex belong to U. Thus, when the union S U U U R 
becomes equal to V(G) for the first time, we have the inequality 2#S < 
#U + ifcR. In particular, it implies that #5 < re/3 at this stage. Let g be 
the girth of G (that is, the length of the shortest cycle in G). When we 
reach the stage S U RU U = V(G), every new soft vertex gives rise to at 
least g — 1 rigid vertices. This means that no more than #U/g vertices will 
be added to S. Therefore, in the case (g > 4), i.e., when G is triangle free, 
the number of soft vertices #5 < n/3 + #E//4 < n/2, as jflJ < 2n/3. Thus, 
the speed of our algoritm is on a par with the fastest algoritms that find a 
single edge coloring or a single Hamiltonian cycle in a cubic graph, or even 
better (cf., e.g., [2]). The absence of short cycles makes the algorithm even 
faster with complexity bounded from above by 2 n ^ 9+2 ^ 39 ~ 2 n//3 . (Note that 
the presence of double edges does not slow down the algorithm because at 
least one of their two common endpoints is rigid.) 

It is instructive to compare the above complexity estimate with the results 
of [3]. Let Ai, . . . , X n £ [—1, 1] denote the eigenvalues of the (normalized) 
adjacency matrix of a simple cubic graph G. Consider the mean// and the 
variance a of the exponents e^' , i = 1, . . . ,n. For each fixed re the points 
(/i, a) form clusters called filars that enjoy a fractal-like structure. From 
the results of [3] combined with the above considerations it follows that the 
closer is the point (/i, a) to the origin, the faster works our algorithm for the 
corresponding graph. 

This algorithm was implemented in C++ code and compiled on a Win- 
dows x86 machine (Pentium IV 3.40 GHz processor with 1 Gb of RAM) 
using CCG GNU Compiler (the program code is given in Appendix). A 
good benchmark for testing programs that search for a Hamiltonian cycle is 
provided by the Horton graph [3], displayed on Fig. 2. 2 It is a cubic bipar- 
tite graph on 96 vertices without Hamiltonian cycles, but with many "long" 
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Figure 2. Horton graph. 



cycles (that is, cycles of length close to 96). Some programs choke when they 
reach such a long cycle, not being able to transform it into a Hamiltonian 
one. Our program completed the search in 5447930319 « 1.268445 x 10 32 
steps (so that the actual complexity is of order 2 n / 3 ). The process took 6336 
sec. of machine time, found 143982592 Tait colorings and no Hamiltonian 
cycles. 

Our program is an open source program and its ANSI C++ code is avail- 
able at the following address: 

http://www.unisanet.unisa.edu.au/staff/homepage.asp?Name= Vladimir. Ejov 



(we do not present it here because of its length). The code does not use 
any platform specific header files, and with minor modifications it can be 
compiled with essentially any C++ compiler that is not mentally challenged. 
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